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Abstract. We investigate the scattering of scalar harmonic source fields by a pe- 
riodic pillar, that is, a spatial structure that is periodic in one dimension and of 
finite extent in the other two. Uniqueness of scattering solutions can be abstracted 
by guided modes. Extending results for periodic slabs to pillars, we give conditions 
under which "inverse" pillars cannot admit guided modes. In addition, we present 
a new construction of guided modes at frequency u and Bloch wavenumber k with 
u embedded in the continuous spectrum for each k. These guided modes have pe- 
riod larger than the period of the structure, and possess a real dispersion relation 
u = oj{n), which is atypical of modes at embedded frequencies. 
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1 Introduction 

Plane-wave scattering and guided modes are important problems in the study of photonic 
crystals and periodic structures. Guided scalar modes are Helmholtz fields exponentially 
trapped within periodic structures in the absence of source fields originating from the ex- 
terior, and this special feature has lead to many applications of photonic crystals such as 
waveguides and light filters. In the case of a periodic slab finite in one direction and periodic 
in the other two directions, variational techniques such as in [1] [5j have been applied to an- 
alyze guided modes. In p[], Bonnet-Bendhia and Starling formulate the scattering problem 
and guided mode problem for periodic slabs and prove the existence and nonexistence of 
guided modes, including guided modes with frequencies embedded in the essential spectrum 
of the corresponding self-adjoint operator. In p], Shipman and Volkov provide a proof of the 
nonexistence of guided modes in inverse structures for photonic crystal slabs, i.e., structures 
that have higher wave speed in the slabs than in the exterior. In this paper, we present a sys- 
tematic mathematical framework for composite structures that are periodic in one direction 
and finite in the other two. We call such structures "periodic pillars" . 

In order to derive the variational formulation, one needs to propose the radiation condi- 
tion through the Dirichlet-to- Neumann operator. In [1], the radiation condition is enforced 
by mapping the Dirichlet boundary data of the Fourier expansion to the Neumann boundary 
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data. Our approach employs Bessel functions (see [7]) to do this, thanks to the fact that, in 
the exterior of pillars, general Helmholtz fields can be expanded as an infinite superposition 
of Fourier harmonics with Bessel functions as coefficients. We propose the radiation condi- 
tion for periodic pillars and formulate the scattering problem through standard variational 
techniques shown in [3] [1] [5] . Naturally, we prove the existence of the solutions to the plane- 
wave scattering problem, and characterize the frequency estimations of guided modes using 
the functional-analytic framework for structures not necessarily piecewise. 

The main original proof in this paper is the existence of guided modes with frequency 
above the cutoff frequency, that is, above the light cone in the first Brillouin zone. The 
existence of guided modes with the frequencies below cutoff are well established in P , and the 
pair of frequency and wavenumber satisfies continuous dispersion relations and an extension 
to periodic pillars is straightforward. The proof of existence of certain embedded guided 
modes is in general difficult. In the proof in section 5.2 of [T], the existence of embedded 
guided modes in periodic slabs is proved. However, this type of embedded guided mode 
can be viewed as non-embedded because the periodicity of the mathematical construction 
is chosen to be larger than the smallest period of the structure and the mode. In fact, 
any non-embedded guided mode in a periodic structure can be treated as embedded if one 
artificially chooses a larger period: this has the effect of reducing the size of Brillouin zone 
in wavenumber space so that the reduced Bloch wavenumber of the mode now lies above 
the light cone. Thus these guided modes are subject to continuous dispersion relations, 
as they are non-embedded from the point of view of their prime period. In this paper, 
we prove the existence of truly nontrivial embedded guided modes in section 4.1. The 
embedded guided modes are obtained by designing a wavenumber-dependent subspace of 
fields that is invariant under the Helmholtz equation and tuning the material parameters so 
that the eigenfunction inside the subspace is located in the regime where all the propagating 
harmonics automatically vanish in this subspace. The eigenfunctions in this subspace are 
automatically guided modes with frequencies embedded in the continuous spectrum for the 
given wavenumber. Our guided modes do not admit a smaller period, yet they persist with 
perturbations of the wavenumbers. The modes admit a continuous dispersion relation above 
the cutoff frequency. 

Our proofs of nonexistence of guided modes are based on the ideas in [T][6]. The first 
proof is an analogy of the proof in |6] for piecewise constant inverse structures. In that study, 
the proof of the nonexistence relies on a restriction on the width of the slabs. The restriction 
on the geometry of the structures is still needed in our proof of nonexistence. We do not 
know whether this restriction is necessary and leave it as an open problem. The second 
nonexistence result is established by using radial monotonicity of material parameters, in 
analogy with Theorem 3.5 of [1] that involves an appropriate Rellich identity. 
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2 Media Structure and Scattering Problem 
2.1 Pillar Structure and Radiation Condition 

Consider an infinitely long "pillar-shaped" structure whose material parameters are periodic 
in the 2;-direction with period 2tt, that is, 

e{x,y,z + 2TT) = e{x,y,z), iJ,{x,y, z + 2tt) = ij,{x,y, z), yx,y,z. (1) 

and bounded in the x, y directions, illustrated by Figure |2.1[ Assume these coefficients to 




Figure 1: a pillar structure periodic in 2;— direction and bounded in x, direction 

be bounded from below and above by positive numbers, e_ < e < e+, /i_ < /i < We 
are interested in time-harmonic fields satisfying the scalar wave equation that are spatially 
pseudo-periodic in z, that is u{x, y, z)e~'^^* , where u satisfies y^z + 27r) = u{x, y, 2;)e^™ 
and K is the Block wavenumber. The field u{x,y,z) is governed by the Helmholtz equation 

V ■ — Vu(a;, y, z) + eoj'^u^x, y, z) = 0, (2) 
/i 

in which e and /x are bounded positive and 27r-periodic in z. 

We use fl = {{x, y, z) : — tt < ^ < tt} to denote one period. Suppose that e = eo, = uq 
ioi r > R and we denote the restricted domain VLr = {{x, y, z) : — tt < z < tt, r = \J x'^ + y"^ < 
R}, whose boundary is = {{x,y,z) G Q : — tt < x,y < TT,r = R} plus the upper and 
lower disks. 



The spatial factor of a time-harmonic acoustic or electromagnetic wave is governed by 
the Helmholtz equation ([2]). By the K-pseudo-periodicity, u can be expanded as an infinite 
superposition of Fourier harmonics in \ Qj^: 

00 00 

u{x,y,z)= J2 E ^-.^W^^''^'^™^'^'' (3) 

m=—oo £=—00 
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where 

{amiHj {rimr) + bmiH^ {Vmr) , if 7^ , 

Cmi + Cm2ln|r|, iir)m = 0,i = 0, (4) 

Cmii\r\^ + Cm£2\r\~^ , if 77^ = 0, £ 7^ 0. 

In this expansion, if wc assume rjm > when 77^ > and irjm > when 77^, the Hankel 
functions {i]mr) are outgoing or exponentially decaying, depending on whether 7]^ is imag- 
inary or real, as r ^ 00, and the Hankel functions Hf{r}mT) are incoming or exponentially 
growing. 

d . . 



The Sommerfeld radiation condition lim r^^^f 



dr 



IK u 



is equivalent to the following 



conditions, which is required for the study of the scattering problem by a periodic pillar. 

Condition 1 (Radiation condition). A field u{r,9, z) satisfies the radiation condition if it 
admits the following Fourier- Bess el representation for r > R: 



meZpUZe l& 



+ E 



I 1/ J 

re 



(5) 



where the subsets Zp^a,e of Z depend on k and are defined by 

m E Zp ^ r]^ > 0,r]m > (propagating harmonics) 
m E Za -r^ rj'^^ = 0,77m = (algebraic harmonics) 
m E Ze 4^ r]^ < 0, —irjm > (evanescent harmonics) . 



i{m+K)z 



2.2 Scattering Problems 

Although this paper deals with guided modes, it would be incomplete without a discussion 
of the scattering problem. In fact, the nonexistence of guided modes is equivalent to the 
unique solvability of the scattering problem. Now consider the specific case of a plane wave. 

Problem 2 (Scattering problem, strong form). Given eo, ei, /xq, A'-i > 0, find u on fl such 
that ^ 

V • —Vu + euP'u = m Q, 
/i 



u is continuous on dfl, 
1 du 



fi dn 



is continuous on dQ, 



u 



E 

m&Zp 



inc i{Kix+K2y+{m+K)z) 



(6) 



u^'^ — u — u^"'^ and its derivatives are K-periodic in z, 
u^'^ — u — u^^'^ satisfies the radiation condition . 
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An incident plane wave e^^'^i^+'^^jz+Ksz) satisfies the Helmhohz equation exterior to the 
pillar, and thus its wave vector satisfies + = rj^. We take k to lie in the first Bril- 
louin zone [—1/2, 1/2). Noting that the function e^^^~^'^ generates the Bessel functions, i.e., 
el = ^^tV^(Z), we let t = e*(''+^«) to obtain e^^^*'^(^+''o) = J^(Z)e*^(^+^»). Then with 
sin 6*0 = —, cos 6*0 = — , and Z = rjmr, the incident wave can be written as a superposition 
of Hankel functions: 

^i{Kix+K2y+K3z) _^i{ri,nr cos sin 6o+rimr sin 9 cos Bo) ^in^z 
_^irimr sm{0+0o) ^i{m+K)z 

= J2MVmr)e^'^'~'''^e^^"^-'-> (7) 
1 

je{0+eo) i{m+K)z 



2 

As a result, the scattering problem of plane waves can be reduced to the linear superposition 
of propagating Fourier harmonics with Hankel functions. 

To analyze the solvability of the scattering problem, it is convenient to reduce to the 
truncated domain flji. Define the pseudo-periodic field space H^^Qr) = {u E H^{Qr) : 

u{x, y, Tc) = u{x, y, — 7r)e^™}. On the boundary F/j, the radiation condition is characterized 

1 1 

by a Dirichlet-to- Neumann map T : {Xr) ~^ Hk ^ (F^j) (as in the Definition 5.19 of [2]) 



where 



i{m+K)z 

ml m,l 




Ime = \ \i\R-\ if m G and £ ^ 



ei 



ii m E Za and £ = 0. 

To satisfy the radiation condition, the harmonics in (|3| with Hj{rimT) for m E ZpVJ Z, 
harmonics {Cmi + Cm2 In |r|)e*^^e**^™"'"'^-'^ for m G Za^t = 0, and harmonics with |r|^ for m G 
Za, £ > all vanish. The radiation condition is hence enforced by 

dnU + Tm = on Tr. (9) 

The operator T has a "nonnegative evanescent" part Te and a "propagating" part Tp-. 

T = Te + Tp. (10) 

^ 1 '^H^{rj!n^R) ^ fmi^ if m G Ze, 

{Tef)^t = < \i\R-^f^g, iimeZa and £ ^ 0, (H) 
I 0, otherwise , 
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(T,/)^, = ^^IIW^^^™^' ''"^^ (12) 
[O, otherwise . 

The variational form of the scattering problem in the truncated domain is 
Problem 3 (Scattering problem, variational form). 



where 



a{u,v) - oo'b{u,v) = f{v),yv e Hlinn) 



a{u,v) = I — Vn ■ Vf H / iTu)v, 



(13) 



b{u,v) = / euv, 
f{v) = -[ [{dnu'^'' + Tu'^'')v]. 

The variational form in problem |3] can be written as 

a{u, v) — uj'^b(u, v) = ci{u, v) + C2{u, v) 



with ci(m, v) = /q^(^ Vm ■ Vf + euv) + ^ f-^^(Tu)v and C2{u, v) = — e(a;^ + 1) J^^ uv. Define 

operators Ci and C2 on HliVLp.) by {CiU,v)Hi(nR) = Ci{u,v) and {C2U,v)Hi{nn) = C2{u,v). 
Because of the coercivity of Ci and the compact embedding of L^{Q) into Hl{Q^), we know 
that the operator Ci is an automorphism and C2 is compact. 

If we denote by Z*"*^ the unique element of HKQr) such that (/*"'^, f )_H'i(nfl) = fi'v), the 
variational form of the scattering problem can be characterized by the following operator 
form 

iCu^v) + iC2U,v) = {n,v),^v G Hli^ln), 

i.e. 

Ciu + C2U = r\ 

The Fredholm alternative theory implies that the nonuniqueness of the solution of this prob- 
lem is equivalent to the singularity of the corresponding homogeneous problem C1U+C2U = 0, 
whose weak form is given by 

a{u, v) - uj%{u, v) = 0, e Hj^inji). (14) 

Theorem 4. The scattering problem has at least one solution, and the set of solutions is at 
most finite dimensional. 
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Proof. From equation ([7]), we can express the incident plane wave as a superposition of 

harmonics ^ - [Hl{r]mr) + H^rj^r)] e^^(^+^o)e^(™+'^)^ with m G Zp. By the Fredholm al- 

e 

ternative, the scattering problem has a solution if and only if 

{f''^,w) = 0, for all w G Null(Ci + C2)\ 

i.e. for all w such that 

a{v,w) -uj%{v,w) = 0,Vt; G 

This w satisfies 

a{w,v) -u%{w,v) = 0,Vf G Hlinji) 

and by the decomposition of T, we know that for all m G Zp, Wm = 0. By the definition 
of Z*"'^, showing = is equivalent to showing that J^^^idn + T)u^"-'^w = 0. This is 

satisfied by the function w above. 

The space of solutions is finite dimensional because Ci is invertible and C2 is compact. □ 



3 Guided Modes 

A guided mode is a solution to the Helmholtz equation in the periodic domain in the absence 



of any source field. In the weak form, it is a solution to the homogeneous equation (14). 
The sesquilinear form 

a'^{u,v)= [ -Vu-Vv+ — f {T'^u)v 
can be split into evanescent and propagating parts, 

a^,{u,v)= [ -Vu-Vv + — [ {T^u)v, 
Jur /^o Jfr 

If the frequency and the wavenumber are assumed to be real, the form (m, v) is Hermitian, 
and we have the following theorem. 

Theorem 5. (Real eigenvalues) If the frequency u is real, then u G H^iytn) solves the 



equation (I4) if and only if 

a'^^iu^v) - a'^{u,v) - u%{u,v) = 0,Vt; G H^inn) (15) 

and if and only if 

'ajiu^v) - uj%{u,v) =0,Wve HliQn), ^^^^ 
{u\rn)m = 0,Vm G Zp. 
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The eigenfrequencies can be obtained by applying the min-max principle to the real form 



in (16). When < y the solutions m of ag(M,f)—(X'^6('U,t>) = 0,Vf G HK^Iji) are guided 



modes since this regime admits no propagating harmonics and hence the second conditions 



in (|16|) are automatically satisfied. When u > ^/j^, to be guided modes, these solutions 



u must satisfy the extra conditions (ulrn)^ = 0,Vm G Zp where Zp is nonempty. We will 
design some periodic structures that admit guided modes in the next section. We have the 
following theorem on properties of the frequencies. The proof is similar to that for periodic 
slabs, for which one may refer to [6] [1] . 

Theorem 6. (Eigenvalues and Frequencies) The problem a'^{u,v) — Xb{u,v) = 0,Vf G 
Hl{Q[i) has a nondecreasing sequence of eigenvalues {Xj}j^i, obtained through the minmax 
principle, 

Xj = sup mf — -, (17) 

dimy=j-i,yci/i(!^fl)"e'^'^\o o[u,u) 

which tend to +oo as j — oo. Moreover, the homogeneous problem a'^{u,v) — uj%{u,v) = 
0,Vf G Hj.{QR) has a nontrivial solution if and only if = Xj{u). This frequency can be 
denoted by uj . 

The problem of scattering and guided modes can be posed equivalently through a positive 
form in the entire strip Q. This formulation leads to the determination of the spectrum of 
the scattering problem, and in particular, the eigenvalues below the continuous spectrum. 
However, the formulation on the truncated domain flu is necessary for the determination of 
the embedded eigenvalues, which will be our focus in the next section. 

We can derive the weak form of the guided modes problem: 

as{u,v) =u%s{u,v),yv e Hl{n) (18) 

where 

as(u,v) = [ -Vu ■ Vv, (19) 



bs{u,v) = / euv. (20) 
Jq 

The associated operator is the unbounded operator 

S^u = --V ■ -Vu. (21) 
e fi 

It is defined on the domain D(S'k) = {m G Hl^{fl) '■ There exists a C such that \as{u,v)\ < 
C^/bs{v, f ), Vf G if^(fi)}. This operator is positive self-adjoint and its eigenvectors and 
eigenvalues are solutions of the guided modes problem. We denote the spectrum of 5'^ as a, 
and its essential spectrum as aess- The following theorem is an adaptation of Theorem 4.1 
of [T] to periodic pillars. 
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Theorem 7. i) a C [-^, +00), where /x+ = sup^ fi, e+ = sup^ e; 




and {Aj(K)} is an increasing sequence 



4 Existence and Nonexistence of Guided Modes 
4.1 Existence 

The focus of this section is to find guided modes with frequency u such that is embedded in 
the continuous spectrum of 5*^. As discussed in the previous section, certain extra conditions 
should be satisfied and hence bring the difiicuhy. 

In [1] , guided modes are proved to exist in a symmetric structure and a periodic slab with 
a finer periodicity. The idea is to consider a closed subspace F on which the operator has 
a cutoff frequency that is greater the cutoff frequency on Hl(Qji), and prove the existence of 
guided modes corresponding to eigenfrequencies lying between these two cutoff frequencies. 
These eigenfunctions are automatically guided modes lying in F because their frequencies 
are below the cutoff frequency, but the frequencies are embedded in the essential spectrum 
of Sk for Hl{flfi). In their proof, the embedded guide modes retain the original pseudo- 
periodicity, but they are simply non-embedded guided modes with a smaller pseudo-period. 
By artificially choosing a larger period, any guided modes with frequencies below the cutoff 
frequency can be seen as embedded guided modes in the same structure with the larger 
period. In this paper, we present a proof of the existence of non-artificial guided modes 
with frequencies embedded in the essential spectrum the operator S'^. We only need the 
parameters e, /i to have smaller period, but the guided modes do not have smaller pseudo- 
period. 

Our newly designed pillar is a periodic structure with period ^ for L > 2 in Z that 
supports guided modes with pseudo-period strictly greater than ^. 

Theorem 8. For any k, in the first Brillouin zone of the structure of period 27i, there exists 
e, /i with period ^ for L > 2 that admits a guided mode with frequencies u lying above the 
cutoff frequency. 



Proof Write u e Hl{nR) Fourier expansion u{r, e,z) = VM„(r,^)e'('"+")^ Given 




M, iV G N with 2M + N + 2 = L, define a nontrivial subspace of Hl{n): 

V = {ue Hl{Q) : Um{r, 9) = 0, if |m - j{2M + N + 2)\<M for some j E Z} (22) 

Therefore, for -M + j{2M + N + 2) < m < M + j(2M + N + 2),the coefficients M„(r, 9) 
are 0, and ioi M + 1+ j{2M + N + 2)<m<M + N + l+ j{2M + N + 2),the coefficients 
Um{r,9) are possibly nonzero. 

We claim that eV C V, fi~^V C V. In fact, let {e)m{r,9) be the Fourier coefficients 
of e. The periodicity of the structure implies that {e)m{r,9) = 0,Vr, 6*, except when m = 




m 
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j{2M + N + 2) for some integer j. For any u eV, ii \m- j{2M + + 2)| < M for some 
j e Z, we calculate the m*'^ Fourier coefficient of eu: 



e 

= y~^(e)j(2M+jV+2)'»m-j(2M+Ar+2) 
j 

— , because Um-j{2M+N+2) — for the field u &V. 



Therefore, eu E V. Similarly, n^^V C V. 

Therefore, the subspace V is also invariant under the operator V ■ ^V. Thanks to the 
invariance properties, wc can consider the Helmholtz equation in the subspace V. The 
solution -u G to the weak formulation a^{u,v) — u'^b{u,v) = 0,\/v E V is also a solution 
to a':^{u,v) — uj%{u^v) = 0,Vf G HK^Ir). In fact, for any field u E V and v G V-^, 
V ■ -Vm + cj'^eu G V implies that V • -Vuv + cu'^euv = for all v G V^. Integrating it we 
obtain 

1_ /"^ 1 /" „ /■ 1 '■^ 



V ■ —Vuv + / e?7,?; = — / dnUV — / — Vit • V^; + / euv 
Jn ^0 Jr^ Jq 

= -a'^{u,v) + h{u,v) 
= 0. 

We can obtain a pair (cu, m) by applying the min-max principle to the Rayleigh quotient 
"^fl^^ on the subspace V to obtain and solving the equation \j{uj) = oJ^ ■ Since u is 

continuous and decreasing from +oo to in ei,/ii seperately, one can choose the material 

parameters such that eo/xoi^^ — + 1 + '^)^ < < eo/xoi^^ — + i-^- P^i^ ('^; 

there are 2M + 1 values — M, — M + 1,...,M — 1,M of m corresponding to propagating 
harmonics. 

The field u obtained in the space V is automatically a guided mode, as the propagating 
harmonics automatically vanish in the subspace V . □ 

As an example, if we let M = iV = 0, then 2M + iV + 2 = 2, 2M+1 = 1, iV+1 = 1, The 
pillar has period tt and e2/+i = for all j, and we can allow one propagating harmonic. Wc 
apply the min-max principle on the space V = {u E Hl{Qji) : U2j = 0,Vj} and by choosing 
proper ei we can obtain an eigenfunction of smallest period 27r that is automatically a guided 
mode. 

If we take M = 1, = 0, then 2M + + 2 = 4, 2M + 1 = 3 and + 1 = 1. Let 
e,n have period 7i/2 and so e.j = for j ^ 4Z, or say Vj, and we can allow to have up 
to 2M + 1 = 3 propagating harmonics. One can minimize the Rayleigh quotient on the 
space V — {u E Hl{ilji) : M4j_i = U4j — 1*4^+1 = 0, Vj}. If we take M — N — 1, then 
2M + A^ + 2 = 5, 2M + 1 = 3, and N + 1^2. The parameters e and /x have period 27r/5 and 
can be allowed to have up to 2M +1 = 3 propagating harmonics. We apply the min-max 
principle on the space V = {u E HI^{Qr) : u^j+i = 1*5^+2 = '^5j+3 — 1*5^+4 — 0,Vj}. The 
pseudo-period of the embedded guided mode is 27t. 
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In our design, the wave number k can be nonzero and there exists a continuous embedded 
dispersion relation The guided mode is robust with respect to k. It is also noticed 

that the modes are subject to the periodicity 2A1+N+2 ■ material is perturbed in a say 

that destroys he small period lose this periodicity while retaining the period 27r, the guided 
mode vanishes. 

This design can also be understood as an existence proof of a guided mode with a larger 
pseudo-periodicity. If we assume the smallest period of the pillar is 27i, embedded guided 
modes with period (2M + + 2)27r can exist. 



4.2 Nonexistence 

Nonexistence results for periodic slabs can be found in [6][T]. In [6], the nonexistence of 
guided modes in inverse structures is discussed. Consider the piecewise constant material 
as in Theorem [6j An inverse structure is a periodic structure with the material parameters 
ei,/ii less than the corresponding parameters eo,/io iii the exterior of the material. The 
proof in [B] requires that the slab satisfy a certain restriction condition. The proof of the 
nonexistence includes introducing the subspace X in which the propagating and linear har- 
monics vanish then estimating the minimum of the Rayleigh quotient. With the restriction 
of the slab width, it is shown that the Rayleigh quotient is strictly bounded below by in 
inverse structures, and hence the weak problem has no solution in X. We use an analogous 
restriction on the radius of the pillar in our proof, and whether this restriction is necessary 
remains an open problem. 

In [1], the assumption is on the parameters only. It is assumed that there exists one 
plane parallel to the slab such that the material parameters e, /x are nondecreasing in the 



direction perpendicular to the slab. In Theorem 10, we present an analogous condition that 
the material parameters be nondecreasing in the radial direction. The proof involves an 
appropriate Rayleigh identity. 

Theorem 9. (Nonexistence of guided modes) Assume that in Q, < e < eo and fi- < fi < 
fio- Let the frequency u and the wave number k, be given in the first Brillouin zone 
Suppose that the radius R of the pillar satisfies 



1 1- 

"2' 2' 



R < ^ (23) 

Then the periodic pillar does not admit any guided modes at the given frequency and wavenum- 
ber. 

Proof. We restrict to the subspace X C HI^Qr) with 

X = {ue H^^inR) : I u{x, y, z)e-'^%-'^"'^^^' = 0, if either m e Zp, or m e and ^ = 0} 



The form a'^{-, ■) is conjugate symmetric in X, and the weak problem (15) is equivalent to 
a'^(u, v) — u^b{u, f ) = on X, as well as a'^{u, v) — uj%{u, v) = Q for all v G X-^ . This gives 

rise to a finite number of extra conditions {dnu\Yif)mi = 0,Vm E Zp 01 m E Z^.i = 0. 
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Consider the eigenvalue problem a^(u,v) — auj%{u,v) = on X. On X, a^(u,v) — 
a'^{u, v). The problem of guided modes is solved by minimizing the quotient fg^ on X. Of 
course, the field u should satisfy the following radiation condition: 



idnu\rn)me + ^rneiu\r 



RJmi 



0, Vm e Ze or m E Za and £ ^ 0. 



(24) 



We first let ei = eo, //i = i^o- The eigenfunctions satisfy a strong form of the Helmholtz 
equation 

(V + iAt)^'0 + aeoHoui'^ijj = in 

Tjjex, rv + 9nV'|r = o (25) 

■0 satisfies pediodic boundary conditions in X. 
In D,R, the separable solutions are in the form of 

A„.J.(|Cn^|r)e^^V(™+'^)^ifC^>0, 



A™,/,(|C™|r)e'^V("^+'^)^ ifC^<0, 
[Cmi + C„2 In \r\] e''^^e^(™+-)^ if (1 = 0, and £ = 0, 
[Cmti\rf + CmMr\-'] e^^^e^(-+«)^ if = and £ 7^ 0, 



(26) 



where — ck€qiiquj'^ — (m + k)^. 

Wc treat the cases for m separately. 

Case 1: m G Zp, i.e. 77^ > 0. In this case, the propagating harmonics should vanish, and 
{u\rji)jne ~ 0- Cm > 0> assume (m > 0, then 

M\Cm\R)^o, 

so ji = QmR = \/ oceoHo — (m + k)'^R, where ji is a zero of Je,{x). The eigenvalues are given 

by 

f + (m + A;)2 



q; 



The Bessel function J; (-2) has a sequence of zeros, and the corresponding a form a sequence 
of eigenvalues {0;^}^^ with all possible ji and m e Z. According to our assumption of the 
radius of the pillar, the eigenvalues 



> 
> 

> 



1 


[1 + 


1 








2 




1 








2 




1 




" 1 




2 





> 1 
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If Cm = 0) the pillar does not support such harmonics for £ = 0. For £ 7^ 0, the separable 
solution Cmiir^ + Cmi2f~^ should satisfy Cmi^R^ + C'm^2-R~^ = which is not possible. If 
Cm < 0, we assume Cm = i\Cm\ and 

h{\UR) = ^. 

It is not possible since the modified Bessel functions Ii have no zeros except at 0. 

Case II: m G Zf., i.e. 77^ < 0, and 77^ = i\rjm\- In this case, the conditions in (24) for m 
should be satisfied. If Cm > 0; and we assume Cm > 0, then 

-^Ji{C,m^)\T=R = —ImlJeiCmR) 

where 7m£ = — ^"^(^''^^j^'' • The value of R can be solved, and by comparing Cm and rj'^, one 
knows that a > 1. If Cm = 0; 'we also have a > 1. If Cm < 0; '^^ assume Cm = ^|Cm|- Then 

\Cm\mUR) + in.M\UR) = o. 

However we know that /^(ICml-R) > 0, ■jme > and Ie{\(rn\R) > 0, and consequently the left 
hand side cannot be 0. 

Case III: rj^ = and £ 7^ 0. The condition {dnu\rji)^f + lme{u\rji)^^ = should be 
satisfied. If Cm > 0; then a > 1. If Cm = 0; there is no solution. If Cm < 0, ICml-^^^dCml-R) + 
ImihiCmR) = 0. It is not possible. 

Case IV: r]"^ = and i = 0. The guided modes satisfy {u\rn)^^ = 0. If Cm > 0; then 
a > 1. If Cm = 0) there is no solution. If Cm < 0, /^(ICIm-R) = 0. It is not possible because 
the Bessel function Jo has no zero. 

In general, when ei = eo, any eigenvalue a > 1. 

In the case when /i_ < /i < /^o and e_ < e < eo, the quotient ^^^^^ with coefficient e, /i 
is greater than or equal to the quotient with quotient eo,/io- Under our assumption of the 
size of the pillar, the number a > 1. As a result, there exists no guided mode because the 
number a > 1 does not correspond to a guided mode. □ 

Theorem 10. Assume there is a pair Xo,?/o such that for all — vr < z < n and any vector 
^0 = ('"ox; '"oj/) 0), the material parameters e,fi are nondecreasing along the direction of Vq, 
that is, the weak directional derivatives Ve ■ Vq, and Vfi-Vo are nonnegative. Then there 
exists no guided mode. 

Proof. Using polar coordinates, we observe that 
du 

V ■ (r—fx-^Vu) =V(rUr) ■ /i" + ruJV ■ fj.~^Vu) 
or 

=rUr{V ■ fx^^Vu) + Ur{Vr ■ jjT^Vu) + rVur ■ jj-^^Vu. 

We integrate this to obtain 

Jvn Jna JflR JflR 

I reUrU + / Ur(Vr ■ /i~^Vn) + / rVur ■ fJ,~^Vu. 
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Adding its complex conjugate, we have 

-1 Dl 



Jtr or Jq^ dr J^^ Jq^ J^^ dr 



We integrate by parts in r for terms including r^, 



er- 



d\u\' 
dr 



2-K /.TT /.i? 




d\u? 
er^r — rdrdzdO 



J -TV JO 

2-K /.TT l-R 




er 



10 J -TV Jo 



dr 
dr 



-drdzdO 



ner^\u\^\^dzde - [ 2er\u\^drdzde - [ r^^\ufdrdzd9 

2e\u\^ - 



■TT 



eR^\u\Xdzde 



r — \u\ 
dr 



and 



[ i_i-^r^ = [ ^ [\-'R''\u\Xdzde - [ 2ii-^\u\^- [ 
Jqh or Jo J-TT Jq.r JOfl 



dr 



The previous identity becomes 



dr 



— — U! 



R\o\u{R)\'^dzd9 - / 2e|VKp 
+ / Ur{Vr ■ ii~Wu) + / Ur{Vr ■ ii~Wu) 



n, 



+ 



27V f'Tl 




R^liQ ^\u{R)\^dzde - I 2ii-^\Vu\ 



J--K 



de, , 
r — \u\ 
dr 



dfT^ 
dr 



\u\ 



Since the field satisfies the Helmholtz equation, we can replace — J^^ fi ^ | Vwp by — a;^ e\u\'^+ 
/Iq^ /r^ wTrM to obtain 



1^0 R\ 



1dI^"|2 



dr 



2iT fn 




J-n 



R'eo\u{R)\''dzde + u' I 2e\Vu\'^ + uo^ f r^\u\ 

Qr Jur or 



JQr JQr 

[ ^ r R''liQ'\u{R)\''dzd0-2uj^ [ e\u\'' + 2ii-' [ uTru - [ 
Jo J-n Jcir Jvr Jqr Or 
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and so 



a /•27r p-iT /•27r p-K 

2 1 /io'^l^l'W / / RhoHR)\^dzde - 



R^fiQ^\u{R)\'^dzde 

[ r^lul"^ + [ u,.(Vr ■ + [ Uriyr ■ ^~^Wu) 
dfi-\ 



|mP + 2fiQ ^ / uTru. 
Jtr 



In this identity, 



du 
dr 



where Vm = + + ^||6/, r = (cos^, sin^, 0), 6> = (- sin^, cos^, 0), 2 = (0,0,1), and 
jVup = + -^l^eP + Simphfy it to obtain 



R 



R 



=2 / fi^^R&' + u 
Fr dr JO 



R 

27V PIT 

2 I , ,2 / / r>2 



R'eo\u{R)\^dzde 



27V PTT 



R'fiQ'\Vu{R)\'dzde 
(27) 



The left-hand side of the identity (27) is nonnegative by our condition on the material 



parameters, and it vanishes if and only if ||M||_H-i(nH) = 0. If we assume m is a guided mode, 
and u has the expansion 



u{r,d,z)= J2 $^a^,ff,i(r^„r)e^^V(™+'=)^+ E 

me2e e meZa ^7^0 



then the terms with m G Za oi the right hand side of (27) is a sum of multiples of 

Since Hj{r]mR) and Hj' {rjmR) are exponentially decaying as i? — )■ 00, in this limit, the 
limit of the right hand side is 0. On the other hand, the left hand side does not converge to 
if M ^ 0. So u = 0. □ 
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